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ABSTRACT 


Lifting-body  configurations  are  constructed  from  stream  surfaces  generated 
by  means  of  supersonic  flows  past  inclined  circular  and  elliptic  cones.  By 
such  means  waverider  shapes  are  devised  with  curved  surfaces  and  known  pressure 
fields  and  shock-layer  structures.  The  conical  flow  fields  stem  from  perturba- 
tions of  the  basic  axisymmetric  cone  flow  arising  from  small  angle  of  attack 
and  small  cross-section  eccentricity.  The  approximate  results  are  analytic 
and  in  the  form  of  hypersonic  small  disturbance  theory.  Various  possibilities 
for  waverider  shapes  are  discussed.  Design  formulas  are  presented  that 
determine  how  the  Mach  number,  angle  of  attack,  cross-section  eccentricity, 
and  characteristic  cone  angle  affect  the  waverider  shape,  pressure  distribution, 
and  shock-layer  structure. 


i 

1 

i 


3 


ACKNOWLEDGEMENTS 


This  Investigation  was  sponsored  by  the  United  States  Air  Force,  under 
Grant  No.  77-3468.  Thanks  go  to  Dr.  Donald  Daniel,  the  project  monitor, 
for  his  Interest  and  encouragement.  Many  enlightening  discussions  with 
Professor  Martin  Jischke  are  also  gratefully  acknowledged. 


TABLE  OF  CONTENTS 


Page 

ACKNOWLEDGEMENTS  Hi 

LIST  OF  ILLUSTRATIONS  vl 

LIST  OF  SYMBOLS  vi  1 i 

1.  INTRODUCTION  1 

2.  FUNDAMENTAL  WEDGE-DERIVED  WAVER IDER  4 

3.  GENERAL  CONSIDERATIONS  FOR  CONICAL  FLOWS  8 

3.1  Coordinates  and  Geometry  8 

3.2  Shock-Layer  Structure  9 

3.3  Conical  Stream  Surfaces  10 

4.  WAVERIDERS  DERIVED  FROM  INCLINED  CIRCULAR  CONES  18 

4.1  Stream  Surfaces  18 

4.2  Approximate  Integration  18 

4.3  Lip  Angle  20 

4.4  Spacing  of  the  Stream  Surfaces  21 

4.5  Waveriders  with  a Freestream  Upper  Surface  22 

4.6  Waveriders  with  a Complementary  Wedge  Surface  25 

5.  WAVERIDERS  DERIVED  FROM  ELLIPTIC  CONES  42 

5.1  Stream  Surfaces  42 

5.2  Approximate  Integration  43 

5.3  Lip  Angle  44 

5.4  Waveriders  with  a Freestream  Upper  Surface  44 

6.  WAVERIDERS  DERIVED  FROM  INCLINED  ELLIPTIC  CONES  50 

6.1  Stream  Surfaces  50 

6.2  Maximum-Entropy  Surface  52 

6.3  Lip  Angle  53 

6.4  Waveriders  with  Freestream  Upper  Surfaces  54 

6.5  Further  Remarks  55 

7.  SOME  CONSIDERATIONS  REGARDING  CONTROL  SURFACES  59 

7.1  Control  Flaps  59 

7.2  Vertical-Fin  Control  Surfaces  59 

7.3  Vertical  Fin  for  the  Elliptic-Cone  Waverider  61 

7.4  Vertical  Fin  for  Inclined-Cone  Waveriders  62 

8.  CONCLUDING  REMARKS  66 


iv 


Page 

REFERENCES  68 
APPENDIX  A:  FORMULAS  FOR  THE  BASIC  AXISYMMETRIC  CONE  FLOW  70 
APPENDIX  B:  FORMULAS  FOR  INCLINED  ELLIPTIC  CONE  FLOW  PERTURBATIONS  72 


V 

«* mmm  ■- jgg* - «»■ 


LIST  OF  ILLUSTRATIONS 


Figure  2.1  a)  Wedge-Shock  Configuration 
b)  Wedge-Derived  Waverider 

Figure  3.1  Cone  Coordinates  and  Geometry 

Figure  3.2  Shock  Eccentricity  Factors,  y = 1.4 

Figure  3.3  Inclined-Cone  Velocity  Perturbations 

a)  radial 

b)  polar 

c)  azimuthal 


Figure  3.4 

Elliptic-Cone  Velocity  Perturbations 

a)  radial 

b)  polar 

c)  azimuthal 

15  | 

1 

I 

!j 

Figure  3.5 

Azimuthal  Velocity  at  the  Body  Surface 

16  :! 

f ID 

Figure  3.6 

Surface  Pressure  Perturbations 

• ”1 

17  ill 

: 2 

Figure  4.1 

Circular-Cone  Stream-Surface  Function 

28  ll 

; n 

Figure  4.2 

Circular-Cone  Lip  Angle 

29  11 

Figure  4.3 

Stream  Surfaces  For  Inclined  Cone 

i 

30  ! 

Figure  4.4 

Circular-Cone  Waverider  With  Freestream  Upper 
Surface:  Positive  Dihedral 

31 

i 

Figure  4.5 

Circular-Cone  Waverider  With  Freestream  Upper 
Surface:  Zero  Dihedral 

32 

Figure  4.6 

Circular-Cone  Waverider  With  Freestream  Upper 
Surface:  Negative  Dihedral 

33 

Figure  4.7 

Circular-Cone  Waverider  With  Wedge-Shock  Upper 
Surface 

34 

Figure  4.8 

a/5  Versus  K.  For  Various  Dihedral  Angles 

0 

3b 

Figure  4.9a 

Wedge  Shock-Circular  Cone  Waveriders, 
a/  5 = 0.1 

36 

Figure  4.9b 

Wedge  Shock-Circular  Cone  Waveriders, 
a/5  = 0.2 

37 

Figure  4.9c 

Wedge  Shock-Circular  Cone  Waveriders, 
r./5  = 0.3 

38 

vi 

— • - -j  - - — 

- ■ ■ ■ ■ ~ 

Page 

1 

Figure  4.9d 

Wedge  Shock-Circular  Cone  Waveriders, 
a/6  - 0.4 

39 

1 

- 

Figure  4.9e 

Wedge/Shock-Circular  Cone  Waveriders, 
a/ 6 * 0.5 

40 

Figure  4.9f 

Wedge  Shock-Circular  Cone  Waveriders, 
a/6  = 0.6 

41 

F 

Figure  5.1 

Stream  Surfaces  for  Elliptic  Cone 

46 

E 

r 

i> 

t 

J, 

Figure  5.2 

Elliptic-Cone  Stream-Surface  Function 

47 

■ 

Figure  5.3 

Elliptic-Cone  Lip  Angle 

48 

‘3 

;-i 

! 

lx 

Figure  5.4 

Elliptic  Cone  Waveriders  With  Freestream 

Upper  Surfaces 

49 

1 4 

1 

jj 

P 

t 

Figure  6.1 

Stream  Surfaces  for  Inclined  Elliptic  Cone 

56 

a 

1 

t 

t 

F 

f. 

t 

Figure  6.2 

Azimuthal  Velocity  Ratios  at  Shock  and  Body 

57 

| 

\ 

Figure  6.3 

Inclined  Elliptic-Cone  Waverider  With 
Freestream  Upper  Surface:  Positive  Dihedral 

58 

? 

l‘ 

i 

Figure  7.1 

Geometry  for  the  Vertical  Fin 

63 

l 

> 

[■ 

k 

Figure  7.2 

Elliptic-Cone  Waverider  with  Fin  Starting 
at  Vertex 

64 

( 

j 

f 

f. 

Figure  7.3 

Elliptic-Cone  Waverider  with  Fin  Starting 
at  Half  Length 

65 

i 

f 

vi  i 

■j  i 

J 

■ TJ 

LIST  OF  SYMBOLS 


r.0,4. 

s 


u„,v 
0 0 


U-|  ,V1  , W-j 


U1  ,v1 ,w1 


a 

& 

Y 

6 

£ 

e, 

p 


tangent  of  semivertex  angle  of  semiminor  axis 

tangent  of  semi  vertex  angle  of  semimajor  axis 

pressure  coefficient  for  basic  circular  cone 

perturbation  pressure  coefficient  for  angle  of  attack 

perturbation  pressure  coefficient  for  cross-section  eccentricity 

eccentricity  factor  for  elliptic  cone 

angle-of-attack  shock  eccentricity  factor 

elliptic-cone  shock  eccentricity  factor 

hypersonic  similarity  parameter, 

freestream  Mach  number 

unit  normal  vector 

spherical  curvilinear  coordinates 

entropy 

radial  and  polar  velocity  components  associated  with  basic  circular 
cone 

spherical  velocity-component  perturbations  associated  with  angle 
of  attack 

spherical  velocity-component  perturbations  associated  with  elliptic 
cone 

angle  of  attack 

unperturbed  shock  angle 

ratio  of  specific  heats 

semi  vertex  angle  of  basic  circular  cone 

perturbation  parameter  for  cross-section  eccentricity 

density  ratio  across  basic  shock 

density 


B/5 


3 


LIFTING  BODIES  DERIVED  FROM 

SUPERSONIC  FLOWS  PAST  INCLINED  CIRCULAR 
AND  ELLIPTIC  CONES 

1.  Introduction 

The  design  of  an  aircraft  for  high  supersonic  flight  that  achieves  desire- 
able  aerodynamic  behavior  and  still  accommodates  such  demands  as  propulsion, 
structures,  materials,  and  operations  is  a very  challenging  task.  A discus- 
sion of  the  requirements  of  hypersonic  aircraft  is  given  by  Kiichemann  [1]  and 
by  Kiichemann  and  Weber  [2].  When  the  Mach  numbers  are  sufficiently  high  that 
the  flow  disturbances  are  Intrinsically  nonlinear,  a treatment  of  the  problem 
by  means  of  linearized  theory  is  not  appropriate.  A generalized  study  of  the 
problem  by  means  of  numerical  solutions  of  the  nonlinear  governing  equations  of 
fluid  mechanics  is  very  formidible  indeed.  For  this  reason  the  few  known  exact 
solutions  for  flows  past  elementary  geometries  are  extremely  Important.  These 
basic  exact  solutions  provide  Insight  and  a knowledge  of  fundamental  physical 
features  associated  with  high-speed  flow.  In  addition,  they  can  also  be  used 
as  building  blocks  for  flows  past  more  complex  geometries. 

A basic  scheme  for  deriving  ex..ct  solutions  for  three-dimensional  lifting 
bodies  by  means  of  simple  two-dimensional  wedge  flows  was  set  forth  by 
Nonweiler  [3],  These  results  were  elaborated  upon  by  Venn  and  Flower  [4],  Nardo 
[5],  and  others.  The  simplest  configurations  thus  derived  are  called  caret 
wings,  or  caret  waveriders,  because  of  their  caret  shape.  Because  these 
aerodynamic  shapes  are  derived  from  basic  two-dimensional  flows,  they  generally 
Involve  flat  surfaces  and  concomitant  sharp  corners  where  these  surfaces  inter- 
sect. These  sharp  corners  may  be  undesireable  when  factors  such  as  viscous  and 
heating  effects  are  taken  into  account. 
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Corresponding  to  the  flat  surfaces  generated  by  the  basic  wedge-shock 
flow,  curved  surfaces  can  be  generated  by  utilization  of  the  stream  surfaces 
associated  with  the  axisymmetric  supersonic  flow  past  a circular  cone.  Such 
surfaces  were  devised  by  Jones  [6]  and  Woods  [7] . These  constructions 
generate  curved  surfaces  and  curved  shocks  that  are  attached  to  sharp  leading 
edges.  The  flows  for  both  the  cone-generated  surfaces  and  the  wedge-generated 
surfaces  can  be  classified  as  conical  flows. 

Besides  the  axi symmetrical  conical  flows,  there  are  also  conical  flows 
associated  with  inclined  circular  cones  and  with  non-circular  cones.  These 
flows  generate  stream  surfaces  that  could  be  attractive  for  constructions  of 
lifting  bodies  with  curved  surfaces.  The  analytical  description  of  these 
flows,  as  contrasted  with  a numerical  description,  generally  involves  a per- 
turbation analysis  of  the  basic  axisymmetric  circular-cone  flow.  A straight- 
forward perturbation  analysis  is  not  uniformly  valid  in  all  the  variables, 
but  fails  in  a vortical  layer  adjacent  to  the  body  surface.  The  pressure 
and  azimuthal  velocity,  however,  are  uniformly  valid  across  the  vortical  layer 
[8,9,10].  This  is  very  Important  because  the  azimuthal  perturbation  velocity 
is  pertinent  in  determining  the  flow  stream  surfaces  and  the  pressure  is 
important  in  determining  the  relevant  forces  on  the  surfaces.  Thus  the  first- 
order  straight-forward  perturbation  expansion,  while  being  suspect  at  first 
glance,  is  thus  pertinent  in  determining  the  stream  surfaces  and  related  forces 
on  waverider  configurations  generated  by  the  perturbation  results.  The  object 
of  this  investigation  is  to  generate  various  waverider  configurations  by  means 
of  arigle-of-attack  and  cross-section  eccentricity  perturbations  of  the  basic 
axisymmetric  cone  flow. 

The  present  study  rests  heavily  on  the  fact  that  approximate  analytical 
expressions  are  available  for  the  perturbed  flows  past  circular  [11,12]  and 
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elliptic  [13,14]  cones  at  small  angle  of  attack.  These  results  allow  for 

an  analytical,  as  opposed  to  a numerical,  investigation  to  be  performed  that 
leads  to  results  that  are  simple  and  easily  understood.  The  effects  of  free- 
stream  Mach  number,  pertinent  cone  angle,  angle  of  attack,  cross-section 
eccentricity,  and  ratio  of  specific  heats  on  the  shock  shape,  shock-layer 
structure,  stream-surface  shape,  and  surface  conditions  can  be  readily  estab- 
lished. Although  the  results  are  approximate,  they  are  accurate  enough  for 
the  parametric  and  design  considerations  that  are  of  primary  concern  here. 

When  the  trends  and  concepts  have  become  clear,  more  precise  and  elaborate 
schemes  of  analysis  can  be  undertaken  for  numerical  accuracy. 

Our  main  concern  In  this  paper  is  the  generation  of  stream  surfaces  that 
can  be  used  as  solid  surfaces  in  lifting-body  waverider  configurations.  We 
shall  outline  how  the  pressure  distributions  can  be  obtained,  but  we  shall  not 
obtain  actual  lift,  drag,  or  moment  results  here.  Such  results,  and  other 
results  of  interest  such  as  reported  by  Squire  [15],  will  be  the  subject  of 
further  research. 
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2.  Fundamental  Uedge-Derived  Waverider 


The  fundamental  wedge-derived  waverider,  sometimes  called  the  caret 
waverider,  is  derived  from  the  basic  exact  flow  field  produced  by  supersonic 
flow  through  a plane  oblique  shock  wave,  as  shown  in  Fig.  2.1a.  The  flow 
is  deflected  by  the  angle  A when  passing  through  the  shock.  The  plane 
shock  wave  is  inclined  at  an  angle  s with  the  freestream  flow.  Let  the 
line  formed  by  the  common  intersection  of  the  planes  parallel  and  perpen- 
dicular to  the  freestream  flow  with  the  plane  of  the  shock  denote  the  leading 
edge  of  a solid  wedge.  The  upper  surface  of  the  wedge  is  parallel  to  the 
freestream  flow  and  the  lower  surface  of  the  wedge,  inclined  downward  at  an 
angle  a with  the  upper  surface,  v the  stream  surface  of  the  flow  that  has 
passed  through  the  shock  at  the  common-line  intersection  denoting  the  leading 
edge  of  the  wedge.  The  flow  field  between  the  lower  surface  of  the  wedge 
and  the  shock  wave  is  parallel  and  uniform  and  has  properties  obtainable 
from  the  oblique-shock  relations. 

The  oblique-shock  relations  are  well  known  and  can  be  expressed  in  ana- 
lytical form  exactly  for  a thermally  and  calorically  perfect  gas.  Owing  to 
the  intrinsic  interaction  of  the  variables,  however,  the  shock  e.igle,  3, 
cannot  be  obtained  explicitly  as  a function  of  the  deflection  angle,  A . 

It  is  thus  useful  to  utilize  the  approximate  results  of  hypersonic  small  dis- 
turbance theory,  valid  for  large  freestream  Mach  numbers,  M^  , arid  small 
deflection  and  shock  angles.  This  shock-angle  relation  originally  obtained 


by  linnell  [16]  is 


sin  e a jz±  . /TrELr  + I 

jin  a 4 / ^ 4 1 (M.  sin  A 


(2.1) 


The  corresponding  relation  for  the  pressure  coefficient  is 
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, . P ~ P~ 

'P  = 1 0 V 2 


CP  = itl  + /fV+1^  ~ 4 ~ 

(sin  a)2  ? / ' 2 ' (M_sin  l)^ 


(2.2) 


In  inviscld  steady  flow,  any  stream  surface  can  be  utilized  as  a solid 
surface  since  no  flow  crosses  a stream  surface.  We  now  wish  to  visualize 
further  stream  surfaces  in  Fig.  2.1a  and  utilize  them  as  solid  surfaces  in 
a new  configuration.  Mark  out  the  axis  AB  aligned  with  the  freestream  flow 
and  lying  in  the  top  surface  of  the  wedge.  The  normal  projection  of  the  axis 
AB  onto  the  lower  surface  of  the  wedge  is  denoted  by  AC.  Any  plane  passing 
through  the  axis  AB  is  aligned  with  the  freestream  flow,  and  any  plane 
passing  through  the  axis  AC  is  aligned  with  the  shock  deflected  flow.  Now 
select  the  points  D and  E on  the  shock  surface,  positioned  symmetrically 
(but  not  necessarily)  on  either  side  of  the  axes  AB  and  AC.  The  two  plane 
surfaces  passing  through  the  pair  of  three  points  ACD  and  ACE  are  stream 
surfaces  in  the  shock  deflected  flow.  The  two  plane  surfaces  passing  through 
the  pairs  of  three  points  ABO  and  ABE  are  stream  surfaces  in  the  freestream 
flow.  If  we  choose  these  four  new  stream  surfaces  as  the  relevant  solid  sur- 
faces of  an  aerodynamic  body,  then  the  new  configuration  appears  as  shown 
in  Fig.  2.1b.  This  aerodynamic  configuration  is  called  a caret  waverider 
since  It  appears  to  ride  on  the  plane  shock  wave  attached  to  its  under  surface. 
The  upper  surface  is  parallel  to  the  free  stream  and  hence  at  freestream 
conditions.  The  under  surface  of  the  caret  shape  is  at  a uniform  pressure 
given  by  formula  (2.2),  and  the  shock  stand-off  angle  is  given  by  formula 
(2.1).  For  a given  caret  waverider  configuration,  the  shock  wave  will  be 


attached  as  shown  In  Fig.  2.1b  only  for  the  particular  Mach  number  determined 
from  Eq.  (2.1)  with  3 and  A given.  The  caret  waverider  is  a lifting  aero- 
dynamic shape  since  the  pressure  is  higher  on  the  under  surface. 
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3.  General  Considerations  for  Conical  Flows 
3.1  Coordinates  and  Geometry 

We  wish  to  determine  the  stream  surfaces  of  conical  flows  past  slender 
elliptic  cones  at  small  angles  of  attack  [13,14].  We  use  spherical  coordinates 
In  a body-fixed  system  shown  In  Fig.  3.1.  Let  a e tan  Qa  and  b = tan  ob  be 
the  semivertex  angles  of  the  semiminor  and  semimajor  axes  of  the  elliptic 
cone.  Then  the  elliptic  cone  Is  described  by 
tan 


tan  0, 


/T+S  cos  Z5 


(3.1) 


where 


tan  Qm  « 


v?  ab 

/FTF 


■ b /1-e  ■ a /T+e 

• vSF  (l-e2)1/4 


. _ b2  - a2 

e 3 ffr+Tr 


(3.2a) 

(3.2b) 


The  parameter  e Is  a measure  of  the  eccentricity  of  the  elliptic  cone.  For 
small  eccentricities  we  can  expand  Eq.  (3.1)  in  a Fourier  series,  the  first 
two  terms  of  which  are 


Gc  * i - e cos  2$  + 0(e2)  , (3.3) 

where 

^ am  + [3  - 2 s1n20m]  sin  29m  + 0(e4)  , (3.4a) 

e = §•  [1  + {15  - 20  s1n20m  + 3 Sln40m  } + 0(eu)]  sin  23m  . (3.4b) 
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The  parameter  c Is  a new  measure  of  the  eccentricity  and  is  the  appropriate 
small  parameter  In  a perturbation  procedure.  The  parameter  6 specifies  the 
semivertex  angle  of  the  basic  circular  cone  about  which  the  perturbation 
expansion  is  performed. 

The  conical  shock  wave  attached  to  the  elliptic  cone  has  the  form  [13,14], 
for  small  a and  e , 


0$  a S + ag  cos  $ - eg  cos  2$  , 


where 


S/5  ■ [ + ] , (3.5) 

6 

and  K,  h Is  the  hypersonic  similarity  parameter.  The  parameters  2f  and  g 
are  the  shock  eccentricity  factors  associated  with  angle  of  attack  and  elliptic 
cone  eccentricity.  They  are  functions  of  and  are  shown  in  Fig.  3.2. 

3.2  The  Shock-Layer  Structure 

Let  u,  v,  and  w denote  the  r,  0,  and  .}>  spherical  components  of  velocity,  and 
let  p denote  the  pressure.  Outside  the  surface  vortical  layer  and  the  viscous 
boundary  layer,  the  variables  have  the  following  expansions  for  small  a and  e: 

u(0,<t>)  a UQ(@)  + aUi(e)  COS  + EUjC©)  cos  2<t> 

v(0 ,<j> ) ■ V ( 0 ) + aVi  (q)  COS  $ + £V  1(0)  cos  2+ 

* (3.6) 

w(o,4>)  = awjlo)  sin  $ + ew^o)  sin  2 <p 

P (© »4> ) 3 Pq(Q)  + up  1 ( 0 ) cos  $ + epi (0)  cos  2* 


The  lowest-order  terms  in  the  expansions,  with  the  subscript  zero,  pertain 
to  the  basic-cone  solution.  The  first-order  terms  with  the  tilde  notation 
pertain  to  the  angle-of-attack  perturbation,  and  the  first-order  terms  without 
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the  tilde  notation  pertain  to  the  elliptic-cone  eccentricity  perturbation. 

The  pressure  and  the  azimuthal  velocity  components,  wj(0)  and  wi(0),  are 
uniformly  valid  across  the  vortical  layer. 

As  calculated  by  Doty  and  Rasmussen  [12],  the  angle-of-attack  velocity 
perturbations  ui/(Vw6),  vt/V^,  and  wi /Va  are  functions  of  0 and  within 
the  fabrication  of  hypersonic  small -disturbance  theory.  Correct  to  the  first 
order  expansion  of  interest,  the  shock-layer  structure  can  be  normalized 
by  the  variable 


a a 


0 - s 
@ - <$ 


(3.7) 


At  the  cone  surface,  & » 0,  and  at  the  shock,  S * 1.  The  angle-of-attack 
velocity  components  are  shown  In  Figs.  3.3a,  b,  c as  a function  of  S and  . 
various  values  of  K..  The  corresponding  elliptic  eccentricity  velocity 
components  [13,14]  are  shown  In  Figs.  3.4a,  b,  c.  The  body-surface  values 
of  the  azimuthal  velocity  components,  w^sj  and  wj (e_ ) , are  shown  In  Fig. 
as  a function  of  Kj.  The  body  surface  perturbation  pressure  coefficients, 
and  cPl»  are  shown  1n  Fig.  3.6  as  a function  of  . 


3.5 


3.3  Conical  Stream  Surfaces 

The  vector  equation  for  a streamline,  V x ds  a 0,  can  be  written  In 

spherical  coordinates  as 

dr  3 rdQ  3 r sin  od» 
uTs.Ty  v(©,4> ) w(0,d) 


(3.8) 


The  conical  stream  surfaces  are  determined  by  the  last  two  terms  of  Eq.  (3.8). 
For  small  angles  and  to  lowest  order  we  have 


do  _ du>  

%T0] 


(3.9) 
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To  lowest  order,  only  the  zeroth-order  polar  velocity,  vQ(e),  for  the  basic 
cone  and  the  uniformly  valid  azimuthal  perturbation  velocities,  w^e)  and  w^e), 
enter  the  analysis.  Even  though  these  velocity  functions  are  known,  it  does 
not  appear  possible  to  obtain  an  integral  of  Eq.  (3.9)  analytically.  In  the 
cases  when  a or  e vanish  separately,  the  variables  in  Eq.  (3.9)  can  be 
separated,  and  Integration  can  then  be  obtained  at  least  by  quadratures. 

The  polar  velocity,  vQ(o),  can  be  approximated  accurately  for  small 
angles  by  [11,12,13,14,17] 

v0(o)  ■ -Vj3(l  - f|)  (3.10) 

Other  formulas  for  the  basic-cone  flow  and  the  perturbed  flow  are  summarized 
in  the  Appendix. 


i. 


eg  cos  24> 


Figure  3.1  Cone  Coordinates  and  Geometry 
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4.  Waveriders  Derived  from  Inclined  Circular  Cones 


4 . 1 Stream  Surfaces 

We  first  consider  the  case  when  e*Q,  that  is,  a circular  cone  at  angle 
of  attack.  When  na0,  we  can  separate  the  variables  in  Eq.  (3.9),  perform 
the  integration  on  3,  and  obtain 


(4.1) 


where 


3*  i s + ag  COS  4>s  , (4.2) 

and  $ is  the  azimuthal  angle  where  the  stream  surface  intersects  the  shock. 

The  function  wi(3)  is  known  [14]  and  shown  in  Fig.  3.3c,  and  thus  the  quad- 
rature in  expression  (4.1)  can  be  evaluated  numerically.  When  3 proceeds 
to  the  cone  surface,  e -*•  5,  the  quadrature  on  the  left  side  of  (4.1)  diverges, 
and  the  left  side  tends  to  minus  infinity.  Correspondingly,  the  angle  4 on 
the  right  side  of  (4.1)  tends  to  zero.  Thus  all  the  stream  surfaces  become 
tangent  to  the  body  surface  at  >=»0,  the  leeward  ray  on  the  cone. 


4.2  Approximate  Integration 

Perusal  of  the  function  w^a)  in  Fig.  3.3c  suggests  that  wi  can  be 
approximated  by  the  relation 


wi (e ) = A + — , 

(4.3a) 

where 

. . -3wi(3)  - wifa) 

M ■ 3-1 

(4.3b) 

B --  - wl(3)] 

(4.3c) 

and  7 = s/i.  When  this  approximation  is  substituted  into  the  integral  in 
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Zq.  (4.1)  and  the  integral  evaluated,  we  obtain 


Ki  ^ C + ^ [ f^l-Oq  + ft2)  In  [|r]-  In  , (4.4) 


where 


fr,  8 _ 1 2.  wi(3) 

M * 2 6V  ’ 

* - f 2^‘(S)  • 


Equation  (4.4)  can  also  be  written  as 

£ 

g*  K1  g* 

r Qj  6 _li  r 3 + f lii 
L a*"-  6'  0 J *•  o*  + 6 ’ 0 1 


tan 


tan  -g- 


tan 


(4.5a) 

(4.5b) 


(4.6) 


Expression  (4.6)  gives  a relation  between  <3/6  and  $ , with  $ , a/6, 
and  K.  as  parameters.  (Recall  that  a s 0/6  is  given  by  Eq.  (3.5).)  The 
parameter  ^ is  always  positive  since  wj (6)  is  negative  when  a is  positive, 
and  vice  versa.  The  dependency  of  fti  on  K6  can  be  determined  with  the  aid 
of  Fig.  3.4.  The  parameter  ft 2 changes  sign  as  varies,  as  might  be  anti- 
cipated from  perusal  of  Fig.  3.3c.  It  can  be  shown  that  [11,12] 


OO 


We  can  thus  write  as 


(4.7) 


2.  q2(o  1 H W|(6|/Va,}  + 2a (j~  - q) 


(4.8a) 


The  dependency  of  k';  on  K,  is  shown  in  Fig.  4.1.  At  K,  * 0 and  K.  =■  0.9, 

w 3 <5  * 9 

the  parameter  ft2  is  zero.  It  is  positive  when  0<K„<0.9and  negative  when 
K ? 0.9. 
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Several  typical  stream  surfaces  (plotted  as  Z tan  3 vs,  ■&,  with  Z 3 constant 
and  tan  9 * 9,  that  Is,  projections  on  a plane  perpendicular  to  the  cone  axis) 
are  shown  In  Fig.  4.2  for  K5  » 1.3,  */&  » 0.2,  and  y ■ 1.4. 

4.3  U_g.  Angle 

We  shall  refer  to  the  portion  of  a stream  surface  adjacent  to  the  shock 
as  the  ljjo  of  the  stream  surface.  The  acute  angle  between  the  lip  of  a stream 
surface  and  the  plane  surface  passing  through  the  lip  and  the  axis  of  the 
body  cone  is  referred  to  as  the  lip  angle,  denoted  by  If  $ - $(s)  denotes 
the  functional  form  for  a stream  surface,  then  the  unit  normal  vector  to  a stream 


surface  Is  given  by  . 

; . 8>  ' slni?3fea 

[1  + S1n28(^)J]‘/2 

where  e and  eQ  are  the  azimuthal  and  polar  spherical  unit  basis  vectors. 
The  lip  angle  Is  determined  by  cos  ^ » e,  • n , and  we  find  that,  for  small 


(4.9) 


angles  9, 


tan  t ■ (&) 


(4.10) 


From  Eq.  (3.9)  we  determine  that 


a* 


aw i ( 3 ) sin  $ 


-t.  , . a s 1 n -!> 
a w i (s)  s 


(4.11) 


Utilizing  Eq.  (4.7)  we  determine  the  Up  angle  to  be 


tan  a ■ f Sin  ^ 


(4.12) 


The  dependency  of  \ on  K{  is  shown  In  Fig.  4,3.  The  ratio  (tan  xJ)csc2  * / (a/6 ) 

Increases  from  zero  to  the  asymptotic  value  5.5  as  K.  increases  from  zero  to 

0 

infinity.  The  lip  angle  thus  becomes  thicker  as  a/-i  and  increase  separately, 
and  It  is  a maximum  when  }.  = 90°. 


1 


r'*1-  **•* 
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4.4  Spacing  of  the  Stream  Surfaces 

In  Eqs.  (4.1)  or  (4.6),  the  value  of  4>s  designates  a particular  stream 
surface,  and  a stream  surface  can  be  represented  functionally  by  $s  = $s(q,$)  * 
constant.  How  closely  the  stream  surfaces  are  spaced  as  Q is  varied  outward 
from  the  cone  surface,  for  a fixed  azimuthal  angle  $,  can  be  evaluated  by  the 
derivative  (3as/30)  with  <|>  held  constant.  From  Eq.  (4.1)  we  obtain 


(ill) 

' 30  1 


. (s)  sin  a 


9 


ov0(a) 


a sin  a 

V 


6 V ©*-<$ 

00 

In  terms  of  the  normalized  shock-layer  variable  S defined  by  expression  (3.7), 


(4.13) 


we  have 


a 

r -T7~ 


sin  a 

&[2+(o-1)&] 


(4.14) 


This  derivative  Increases  when  x increases  and  goes  to  Infinity  when  4 goes 
to  zero  (0  -*•  5),  and  the  stream  surfaces  become  more  closely  spaced  correspond 
Ingly.  As  seen  from  Fig.  3.3c,  the  derivative  Is  Insensitive  to  variations 
In  when  Is  about  0.7.  When  Sis  less  than  0.7,  the  derivative  decreases 

as  Increases,  and  vice  versa  when  % is  greater  than  0.7,  but  to  a lesser 

degree. 


•X— 
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4 . 5 Waveriders  with  a_  Freestream  Upper  Surface 

Any  stream  surface  just  described  can  be  utilized  as  a solid  surface, 
but  a complementary  surface  remains  to  be  described  in  order  to  fashion  a 
closed  aerodynamic  body.  In  this  section  we  choose  such  surfaces  that  are 
parallel  to  the  free  stream.  We  first  demarcate  the  axis  passing  through 
the  vertex  that  is  aligned  with  the  free  stream.  This  axis  Is  inclined 
by  an  angle  a with  the  cone  axis.  Any  plane  passing  through  this  freestream 
axis  is  parallel  to  the  free  stream.  We  select  pairs  of  these  freestream 
planes  that  pass  through  thelipshcck  intersection  of  the  conical  stream 
surfaces,  and  a closed  aerodynamic  body  is  thus  formed.  The  upper  surfaces 
are  pairs  of  freestream  planes  passing  through  the  freestream  vertex  axis, 
and  the  lower  surface  is  a stream  surface  of  the  inclined  circular-cone 
flow  field.  These  surfaces  Join  together  at  the  lip-shock  intersection. 

An  example  is  shown  in  Fig.  4.4  for  which  the  lip-shock  intersection 
occurs  at  $ - 90°. 

There  are  an  Infinite  number  of  such  aerodynamic  waveriders,  depending 
on  how  the  freestream  planes  are  selected.  As  a step  towards  distinguishing 

between  the  different  configurations,  let  us  determine  the  lip  angle,  X , 
between  the  freestream  planes  and  the  conical  stream  surfaces.  The  normal 
unit  vector  to  a freestream  plane  intersecting  the  shock  at  ♦ is,  for  small 
angles 

Co  - (f)O-g)  cos  *s]  e.j  - (j)  sin  e0 


/[o-  (£)(!-$)  cos  4S]2  + (£)2  sin2't>s 


(4.15) 
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The  lip  angle,  £ , Is  determined  by  the  condition 


cos  x ■ n 
w ■ 


(4.16) 


where  n is  the  unit  normal  vector  of  the  conical  stream  surface  evaluated  at 
the  shock,  given  by  Eq.  (4.9).  We  find  for  small  ratios  a/o  <<  1 that 


tan  ’ f c 


1 


o - j (1-g)  cos  $s 


] sin 


(4.17) 


If  the  waverider  lip  angle  is  ever  to  be  a cusp,  then  we  must  have  X * 0. 

w 

The  right  side  of  Eq.  (4.17)  can  vanish  only  when  sin  $ ■ 0 or  when 

f cos  *s  “ ‘ (4-18) 

The  right-hand  side  of  Eq.  (4.18)  vanishes  when  - 0 and  Increases  to  0.87 
when  * od.  When  ■ 0.5  and  a/5  - 0.5,  then  Eq.  (4.18)  yields  $ ■ 0, 
and  a cusp  can  occur  only  at  $s  « 0.  When  is  larger  than  0.5  and  a/5  Is 
less  than  0.5  a cusp  cannot  occur  except  at  $s  * 0°  or  180°.  These  are  generally 
the  conditions  of  interest  and,  as  such,  the  possibility  of  a cusp  is  of  no 
concern. 

From  structural  or  heating  considerations,  the  condition  where  the  lip 
angle  £w  is  a maximum  is  of  interest.  Setting  the  derivative  of  xw  with 
respect  to  -js  equal  to  zero  yields,  for  small  a/5, 


cos  $s 


siil  i 

a 5 


(4.19) 
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When  the  right-hand  side  Is  less  than  minus  unity,  no  relative  maximum  occurs. 

When  the  right-hand  side  Is  small,  we  obtain 

, (*.20) 

which  Illustrates  that  the  maximum  value  of  occurs  when  $s  > it/2.  When 

K5  ■ 1.3,  y “ 1-4,  and  a/4  ■ 0.2,  ■ 8.6°  Is  a maximum  when  $s  a 96.8°. 

This  particular  waverlder  Is  similar  In  shape  to  that  shown  In  Fig.  4.4,  for 
which  $ ■ 90a  and  - 3.57°.  The  upper  surfaces  In  these  cases  appear  to 
have  a positive  dihedral  angle. 

Another  waverlder  can  be  formed  in  which  the  two  upper  freestream  planes 

are  parallel,  that  Is,  the  upper  surface  Is  flat.  An  example  is  shown  In  Fig.  4.5. 

From  geometrical  considerations,  this  situation  occurs  when,  for  small  a/s, 

cos  *s  ■ . (4.21 ) 

For  K3  ■ 1.3,  v “ 1.4,  and  at/4  ■ 0.2,  we  get  $s  ■ 81.4°  and  ■ 8.2°.  This 
waverlder  Is  akin  to  the  waverlder  formed  from  a half  cone  at  zero  angle  of 
attack  with  a symmetry  plane  through  the  cone  axis  Identified  as  a flat,  zero- 
thickness delta  wing.  The  waverlder  shown  in  Fig.  4.5,  however,  while  having 

a flat  upper  surface,  has  a faired  under  surface  with  a "wing"  of  finite 
thickness.  This  waverlder  can  be  said  to  have  zero  dihedral  angle. 

Other  waveriders  in  this  family  can  be  formed  that  have  negative  dihedral 
angles.  These  exist  when  $ is  less  than  the  value  given  by  Eq.  (4.21).  An 

example  Is  shown  In  Fig,  4.6  for  $ ■ 70*  and  x'  ■ 7.6°.  All  of  the  waveriders 

In  this  section  are  lifting  bodies  since  the  lower  conical  surface  is  at  a 
higher  pressure  than  the  upper  freestream  surfaces. 
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4.6  Waveriders  with  a Complementary  Wedge  Surface 

Another  means  of  deriving  a closed  aerodynamic  shape  in  conjunction  with  a 
conical  stream  surface  is  to  use  plane  surfaces  passing  through  the  cone  axis. 

In  these  cases  the  freestream  flow  must  be  deflected  by  the  angle  a.  We  can 
thus  use  the  formulas  for  the  fundamental  wedge-derived  waverider  obtained 

in  Section  1,  setting  a ■ a.  Let  the  angle  between  the  two  cone-axis  plane 

surfaces  be  denoted  by  2< \>,  where  \\i  Is  the  dihedral  angle,  as  shown  In  Fig. 

4.7.  The  plane  (or  wedge)  shock  wave  across  the  top  of  the  waverider  Is 
oriented  at  an  angle  sw  with  the  free  stream,  and  It  Is  related  to  the 
cone  shock  by  the  relation 

aw  - a ■ 0*  cos  ij; 

(4.22) 

■ (s  + COS  $s)  COS  '() 

We  now  note  that  j>s  ■ it  - ip  and  that  aw  Is  given  by  Eq,  (2.1),  with  sin  3 - 3W 

and  sin  a •*  a.  We  thus  rewrite  Eq.  (4.24)  as 


[|  - j g cos  1/]  cos  v 


(4.23) 


Recalling  that  a = 3/3  Is  given  by  Eq.  (3.5),  we  note  that  Eq.  (4.25)  provides 

a relation  between  a/3,  K^,  and  >ii.  Solving  for  a/3,  we  obtain 

0 


a 

3 


Ba  cos  ili  ± J (-^j—  a cos  ) 4 - ( (-'-jp-)  ' - B*- } 

(^)2  - B2 


(4,24) 


where 


9 3 ^-4-  - ^ COS  i Ij 


(4.25) 
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Figure  4.8  shows  a/6  plotted  as  a function  of  for  various  values  of  the 
dihedral  angle  p and  for  y ■ 1.4.  We  observe  that  a/4  Is  double  valued  for 
fixed  values  of  and  p,  but  this  Is  only  apparent  for  u > 75°  since  the 

curves  are  plotted  for  only  the  realistic  values  of  a/5  < 1.  The  condition 

a/6  ■ 0 corresponds  to  zero  flow  deflection  over  the  upper  surface;  the 
dihedral  angle  p Is  such  that  a plane  Mach  surface  sits  on  the  lip  of  the 
waverlder,  and  the  Ups  extend  to  the  cone  surface  as  a pair  of  Infinitesimally 
thin  delta  wings. 

The  lift  on  the  waverlders  In  this  section  may  be  positive  or  negative 

depending  on  the  value  of  a/4.  The  pressure  In  the  conical  flow  field  Is 

given  by 

Cp  (9,+)  Cp  (®)  £p  (0) 

— -J2 • — Jr-  +f— \ cos  p + 0(a2)  , (4.26) 


and  the  pressure  In  the  wedge  flow  field  Is  given  by  (from  Eq.  (2.2)) 


2 

+ 


(4.27) 


The  minimum  pressure  on  the  conical  stream  surface  occurs  at  the  symmetry  plane, 
0 " 4 and  p ■ 0.  The  minimum  lifting  pressure  differential  Is  thus  given  by 


CD  («,0) 

Hc 

T2 


We  also  note  that  [ , , ] 

CB  (4) 

Po  . i , a2  In  a2 

i7  ' 1 + a2  - T" 


(4.28) 


(4.29) 


27 


Further,  from  Fig.  3.6,  we  note  the  C (5)/5  * -4.  A lower'  bound  for  the  value 

Pi 

of  a/4  corresponding  to  zero  lift  can  now  be  found  by  setting  AC  ■ 0 : 

pm 


a2  In  a2 


(4.30) 


This  relation  provides  a/5  for  zero  lift  as  a function  of  and  Is  shown  In 
Fig.  4.7.  The  waverlders  described  by  conditions  below  this  line  are  lifting 
when  the  conical  surface  Is  underneath,  and  the  waverlders  described  by  condi- 
tions above  the  line  are  lifting  when  the  wedge  surface  Is  underneath. 
Approximately,  the  zero  lift  conditions  occur  when  a/4  » 0.4  for  K6  > 1 . 

A systematic  variation  of  the  cone-wedge  waverlder  cross  section  geometries 
Is  shown  In  Figs.  4.9a,  b,  c,  d,  e,  f for  a/5  ■ 0.1,  0.2,  0.3,  0.4,  0.5,  0.6. 

For  each  value  of  a/5,  dihedral  angles  of  1/  ■ 45* , 60°,  70°,  75®,  80°  are 
shown  along-wlth  their  corresponding  values  of  as  determined  from  Eqs.  (4.25) 
or  (4.26).  For  a/4  ■ Q . 1 , 0.2,  and  0.3,  the  bodies  are  lifting  when  the  conical 
surface  Is  underneath.  At  a/5  « 0.4,  the  body  Is  nearly  at  zero  lift.  For 
a/4  ■ 0.5  and  0.6,  the  bodies  are  lifting  when  the  wedge  surface  Is  underneath. 
The  shock  lies  closer  to  the  body  when  Is  larger,  and  hence  when  -p  Is 

larger.  The  standard  conditions  K ■ 1.3,  •<  ■ 1.4,  and  a/5  ■ 0.2  are  represented 
In  Fig.  4.7,  In  which  case  'p  * 57°  and  a ■15’. 


Circular-Cone  Waverider  With  Freestream  Upper 
Surface:  Positive  Dihedral 
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4.9b  Wedge  Shock-Circular  Cone  Waveriders 
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5.  Waverlders  Derived  from  Elliptic  Cones 
5 . 1 Stream  Surfaces 

We  now  consider  the  case  when  a ■ 0,  that  Is,  an  elliptic  cone  at  zero 
angle  of  attack.  When  a ■ 0,  we  can  separate  the  variables  In  Eq.  (3.9), 
perform  the  Integration  on  s,  and  obtain 


where 

<3*  a 3 - eg  cos  s , (5.2) 

and  where  we  have  replaced  6 that  appears  In  Eq.  (3,10)  by  0 , the  elliptic 

c 

cone  angle, 

s 3 ~ e cos  24  . (5.3) 

The  variation  of  4 that  occurs  in  0C  Is  to  be  ignored  In  the  Integration. 

These  operations  are  consistent  In  our  present  first-order  accuracy.  The  azi- 
muthal velocity  wi(0)  Is  known  [13,14]  and  Is  Illustrated  In  Fig.  3.4c.  The 
quadrature  thus  can  be  evaluated  numerically.  When  0 proceeds  to  the  cone 
surface,  Q ■*  0 , the  quadrature  on  the  left  side  of  Eq.  (5.1)  diverges,  and  the 
right  side  tends  to  minus  Infinity.  Correspondingly,  when  the  shock  location 
of  the  stream  surface,  4S.  Is  less  than  90°.  the  angle  4 tends  to  zero.  When 
4S  Is  greater  than  90°,  the  angle  4 tends  to  180°.  Thus  the  stream  surfaces 
that  begin  at  the  shock  for  4$  < 90°  become  tangent  to  the  cone  surface  at 
4 a 0,  and  the  stream  surfaces  that  begin  at  the  shock  for  4S  > 90°  become 
tangent  to  the  cone  at  4 ■ 180°.  The  stream  surface  for  4 * 90°  Is  a 
symmetry  plane  that  Is  perpendicular  to  the  body  at  4 * 90°.  Examples  of 
stream  surfaces  are  shown  In  Fig.  5.1. 
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5 . 2 Approximate  Integration 

Figure  3.4c  suggests  that  w^e)  can  be  approximated  by  the  relation 


Qr  B 

wi(e)  ■ A + -±r- 


where 


A = 


a W^O  ) - w:(9  ) 


-v 

o 


T 


C«i(ec)  - «i(Qs)] 

B s — yn 


(5.4a) 

(5.4b) 

(5.4c) 


where  o = Qs/oc  ■ o + 0(e),  0S  * s + 0(c),  and  ec  ■ 5 + 0(e).  The  errors  of 
order  e In  Eqs.  (5.4)  will  subsequently  be  neglected.  When  this  approximation 
Is  used  to  evaluate  the  integral  In  Eq.  (5.1),  we  obtain 

0 - 0„  0+0, 


ki  ferirh  + ka  In  k»rr]  - (ki  + k2>  1"  C|r]  ■ 'n  , (5.5) 

''c  c r c + c 


where 


k,  • 

u 1 w 

k2  = - f wi(p). 

Equation  (5.5)  can  also  be  written  as 


(5.6a) 

(5.6b) 


- 0 


rt*  *1  0 + 0 


■]  c 


0*  ^2 


tan 


e*  + 0 


tan 


(5.7) 


Expression  (5.7)  Is  a relation  between  ©/i  and  <*  , with  $ , e/s  , and 
as  parameters.  The  parameter  Is  always  positive  since  WjU)  Is  negative 
when  c is  positive,  and  vice  versa.  Figure  3.5  Illustrates  the  dependency  of 
k,  on  K . The  parameter  k2  changes  sign  as  K varies,  which  is  suggested  by 

Q 0 

the  behavior  of  Wj(e)  shown  In  Fig.  3.4c.  It  can  be  shown  that  [13,14] 

3 - 2g/a2  . (5.8) 
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We  can  now  write  k2  as 

4g  + 0(0+1 ) {Wi(fi)/V  } 

k^-I ^T) “ (5 

Figure  5.2  shows  how  k2  depends  on  . The  parameter  k2  vanishes  at  « 1.9. 
In  the  range  0 £ < 1 .9,  k2  is  positive,  and  it  is  negative  when  K$  ^1.9. 


5.3  Lip  Angle 

The  Up  angle  Is  defined  by  cos  x 
surfaces  is  determined  by 


tan  x 


...  ew,(fi)°sin2*s 




e 2g  sin  2«j«s 
6 o(oz-1  ) 


e$*n,  and  for  the  elliptic  conical 


(5.10) 


Figure  5.3  shows  the  dependency  of  x ori  Kfi.  The  ratio  (tan  x)(csc  2<|>s)/(e/6) 
Increases  from  zero  to  the  asymptotic  value  8.7  as  K.  Increases  from  zero  to 
Infinity.  The  lip  angle  becomes  thicker  as  and  e/6  increase  separately. 
The  lip  angle  Is  a maximum  when  <ts  “ 45°  and  1 35° , and  It  is  zero  at  the 
symmetry  planes  $s  = 0°  and  90°. 


5 . 4 Waverlders  with  Freestream  Upper  Surface 

The  axis  of  the  elliptic  cone  at  zero  angle  of  attack  Is  parallel  to  the 
free  stream,  and  any  plane  that  passes  through  the  cone  axis  Is  parallel  to  the 
free  stream.  We  choose  pairs  of  plane  surfaces  that  pass  through  the  cone 
axis  and  Intersect  the  shock  at  the  lips  of  a conical  stream  surface.  When 
the  plane  surfaces  are  the  upper  surfaces  and  the  conical  surface  Is  the  lower 
surface,  a lifting  aerodynamic  waverlder  Is  formed.  Two  examples  are  shown  In 
Fig.  5.4  for  K5  - 1.3,  c/6  = 0.1,  and  X = 1.4.  These  configurations  correspond 
to  lip  positions  of  f = 100°  and  110°.  The  respective  lip  angles  are  found 
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to  be  | x | = 2.2°  and  4.0°.  These  elliptic-cone  shapes  are  similar  in  form  to 
the  circular  cone  shape  shown  in  Fig.  4.6.  The  elliptic  cone  shapes  are  flatter 
on  the  bottom  and  have  thinner  lip  angles. 

The  pressure  coefficient  on  the  undersurface  is  given  by 


Cp(©,*) 


c„  (0) 


~P~ 


+ <T> 


v9) 


cos  2$ 


(5.11) 


When  Eq.  (5.7)  Is  used  to  determine  Q as  a function  of  4>,  the  surface  pressure 
as  a function  of  $ can  be  determined. 
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6.  Waveriders  Derived  From  Inclined  Elliptic  Cones 


6 . 1 Stream  Surfaces 

The  general  stream-surface  structure  for  flow  past  an  inclined  elliptic 
cone  can  be  studied  by  rewriting  Eq.  (3.9)  In  the  form 


where 


(1  -"-Wa  +T4u7 


u = cos 


(6.1) 

(6.2a) 


ka  5aWi(0)1+"lewi(O)  ’ (6.2b) 

k4  +^2cwi(0)  • (6.2c) 

The  difficulty  with  the  Integration  of  Eq.  (6.1)  comes  from  the  combination 


w*(0,u)  s k3  + k4u 


(6.3) 


which  may  vanish  In  the  shock  layer.  The  azimuthal  velocity  vanishes  between 
the  shock  and  the  body  when  w*  * 0.  The  azimuthal  location  $o  where  w*  vanishes 
Is  determined  by 


cos  <j>0  3 


* 

w> 


(6.4) 


The  surface  for  w*  3 0 (dashed  curve)  is  shown  In  Fig.  6.1  for  3 1.3, 

Y a 1.4,  a/$  3 0.2,  and  s/<5  3 0.1.  The  position  where  this  surface  Inter- 
sects the  body  surface  Is  denoted  by  3 *oc  and  is  determined  by 

a *i(0c> 

cos  v ■ - 27  wn^y 


For  the  conditions  In  Fig.  6.1  we  have  $ 3 127°.  The  surface  for  w*  3 0 

lies  on  the  windward  side  of  the  body.  Between  the  body  and  the  surface  w*  3 0 
the  azimuthal  velocity  is  positive,  and  It  Is  negative  between  the  w*  = 0 surface 
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and  the  shock.  For  other  conditions,  cos  $QC  can  be  calculated  with  the 
use  of  Fig.  6.2,  which  shows  the  ratio  w1(5)/wi(5)  as;  a function  of 
for  y ■ 1.4.  Also  shown  is  the  ratio  at  the  shock,  Wi ( 3) /wi (s ) . For  a given 
value  of  K^,  the  ratio  w1(e)/iwi(Q)  Increases  from  the  shock  value  to  the 
surface  value. 

We  can  obtain  a straightforward  approximate  Integration  of  Eq.  (6.1) 
by  treating  ka  and  k4  as  constants.  Experience  shows  that  the  best  values 
to  use  for  these  constants  are  the  values  at  the  body.  For  this  choice, 
the  body  surface  conditions  are  correctly  given,  which  Is  especially  important 
for  the  position  of  the  Intersection  of  the  windward-leeward  dividing  stream 
surface  with  the  body,  denoted  by  $ . 

With  the  left  side  of  Eq.  (6.1)  evaluated  as  for  the  previous  cases, 
we  obtain  the  approximate  Integral  of  Eq.  (6.1)  to  be 


0 - 0„ 

0*  ' 

ki+ki 

“0  + 0 

0*  ’ 

C 

s 

c 

s 

t 

1 

O 

0 

g*  + 0„ 

0 ’ 

L 5 c 

s C 

k5  k6 


"tan 

Ri  esc  $ + ki  cot  4> 

*s 

tan  f 

Ri  esc  'p$  + kx  cot  $s 

where  k5  = ^/(kj  - kx)  , (6.7a) 

k6  = k|/ ( kx  - kx)  , (6.7b) 

0C  = 5 - e COS  2$  , (6.7c) 

0*  * a + ag  COS  4>s  - eg  COS  2$s  . (6 . 7d ) 


Equation  (6.6)  reduces  to  Eq,  (4.6)  or  (5.7)  In  the  separate  limits  e - 0 
and  a ■+  0.  If  the  conditions  are  such  that  kx  = kx,  then  the  appropriate 
limit  must  be  taken  on  the  right-hand  side  of  F.q.  (6.6),  which  corresponds 
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to  setting  k3  * k4  In  the  original  differential  equation  (6.1). 

Figure  6.1  shows  stream  surfaces  generated  by  Eq.  (6.6)  for  the  conditions 
K5  ■ 1.3,  y * 1.4,  a/5  ■ 0.2,  and  e/5  ■ 0.1.  The  dividing  stream  surface 
which  separates  the  flow  that  proceeds  toward  the  leeward  ray  from  the  flow 
that  proceeds  toward  the  windward  ray  occurs  at  <j>  ■ 127°.  This  dividing  sur- 
face Is  a plane  surface  In  this  approximation.  The  actual  stream  surfaces 
should  appear  slightly  different  near  the  dividing  stream  surface  and  toward 
the  windward  ray.  In  this  region  the  actual  stream  surfaces  have  a negative 
value  of  azimuthal  velocity  at  the  shock  and  hence  should  slope  towards  the 
leeward  direction  before  curving  Into  the  radial  direction  at  the  surface 
w*  - 0 and  then  curving  further  towards  the  windward  ray  as  shown.  The 

slope  (3$/d0)  should  actually  be  positive  In  the  windward  region  whereas 
*s 

the  approximation  (6.6)  yields  a slightly  negative  slope.  Because  the  details 
of  the  azimuthal  velocity  field  have  not  been  taken  Into  account,  the  dividing 
stream  surface  and  Its  Intersection  with  the  body  surface  are  not  precisely 
described  by  the  plane  surface  shown,  even  though  the  location  of  the  body 
Intersection  Is  correct.  Both  the  actual  slopes  and  the  approximation 
slopes  at  the  shock  are  small  In  the  windward  region,  and  hence  the  descrip- 
tion shown  In  Fig.  6.1  Is  qualitatively  valid  except  for  the  aforementioned 
discrepancies . 

6.2  Maximum-Entropy  Surface 

Outside  the  vortical  and  viscous  boundary  layers,  the  entropy  has  the 
expansion 

s(Q,'j>)  n SQ ( Q ) + as  i cos  $ + esi  cos  2cf>  , (6.8) 

where  Si  and  s ; are  constants  associated  with  the  angle-of-attack  and 
eccentricity  perturbations.  Let 

S?  = aSj  cos  ^ + e S i cos  2sj>s  (6.9) 


at-  ir 
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denote  the  entropy  perturbation  at  the  shock.  It  can  be  shown  [11,13]  that 
if  2 V ^ 

■ -y(y-1  ) 3 (1-C0)  JZJJJ  [a(l-§)  cos  4>s  + Eg  cos  2*s]  , (6.10) 

where  5Q  ■ p„/p0(3)  and  (o ) ■ yP0(3)/pq(3) . The  maximum  entropy  pertur- 
bation at  the  shock  occurs  where  the  derivative  with  respect  to  «s  of  Eq.  (6.10) 
vanishes.  This  occurs  where 


cos  <p 


sm 


4 eg 


(6.11) 


If  the  value  of  the  right  side  Is  less  than  or  equal  to  minus  unity,  the  maximum 

occurs  at  the  windward  ray,  4>s  • 180°.  For  the  conditions  of  Fig.  6.1 

(K,  ■ 1.3,  y ■ 1.4,  a/5  ■ 0.2,  c/5  ■ 0.1),  the  maximum-entropy  stream  surface 
originates  at  the  shock  at  $s  "133°.  This  Is  on  the  windward  side  of  the 
dividing  stream  surface  which  Is  located  In  Fig.  6.1  at  $s  ■ 127°.  Thus  the 

maximum-entropy  stream  surface  does  not  wet  the  body  surface  but  lies  In  the 

windward  part  of  the  shock  layer.  That  the  body  surface  and  maximum-entropy 
surface  are  not  necessarily  Identical  Is  also  true  for  hypersonic  blunt-body 

flows  (see  Hayes  and  Probsteln  [18],  page  399).  Melnik  [10],  on  the  other  hand, 
In  his  analysis  of  conical  flows  assumed  that  the  body  was  a maximum  entropy  sur- 


face . 


fhe  conical  result  holds  for  all  values  of  K, 


The  ratio  of  cos  a 


sm 


and  cos  d , given  by  Eq.  (6.5)  Is 


cos  * 


sm 


COS  $ 


oc 


Tg 


W i ( 5 ) 


(6.12) 


For  y * 1.4  this  ratio  Is  greater  than  unity  for  K.  > 0,  and  hence  we 

O 


conclude  that  * > a for  K. 

ITl  Q C g 


0. 


6 . 3 U£  Angle 

Corresponding  to  Eqs.  (4.12)  and  (5.10),  we  determine  the  lip  angle, 
for  the  shock-stream  surface  intersection  for  an  inclined  elliptic  cone 
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to  be  given  by 


tan  x ■ tan  x + tan  x 

,2  " 


af  sfSn  s1n  *s  + f n^rrs1n  2*s 

When  the  angles  are  small,  the  lip  angle  for  an  Inclined  elliptic-cone  stream 
surface  Is  equal  to  the  sum  of  the  lip  angles  for  the  corresponding  inclined 
circular-cone  stream  surface  and  the  aligned  elliptic-cone  stream  surface. 


(6.13) 


6 . 4 Waveriders  With  Freestream  Upper  Surfaces 

As  before,  waverider  configurations  can  be  formed  by  utilizing  freestream 
plane  surfaces  that  pass  through  the  vertex  along  an  axis  inclined  at  an  angle 
a with  the  cone  axis,  and  Intersect  the  shock  at  the  lips  of  an  underneath 
conical  stream  surface.  As  an  example  we  choose  the  conditions  shown  In 
Fig.  6.1  with  the  lip  located  at  *s  ■ 90a.  The  waverider  so  formed  is  shown 
in  Fig.  6.3.  The  freestream  surface  Intersecting  the  shock  at  4>s  has  a unit 
normal  vector  given  by 

Co  - (f)(l-g)  COS  * - (f)  g cos  2*1  4 - | sin  * e 

n-  - j =1 Z t==  - = 9 • (6-14) 

/Co  - (|)(l-g)  cos  *s  - (|)  g cos  2*$]2+  (|)2  s1nz*s 

A A 

The  waverider  lip  angle,  x.„,  Is  determined  by  the  condition  cos  xw  ■ n • n, 
where  n is  given  by  Eq.  (4.9).  For  small  a/5,  we  obtain 


tan  X »<  tan  x 2 2 

w o - (|)(l-g)  cos  *s  - (|)  g cos  2*s 

The  waverider  lip  angle  In  Fig.  6.3  is  found  to  be  x * 8,9°.  This  waverider 

W 

has  a positive  dihedral  and  is  analogous  to  the  circular  cone  waverider 
shown  in  Fig.  4.4.  The  elliptic-cone  waverider  in  Fig.  6,3  is  more  shallow 


(6.15) 


55 


and  flatter  on  the  bottom  than  Its  circular-cone  counterpart  In  Fig.  4.4. 

It  also  has  a slightly  greater  span  since  the  attached  shock  Is  broader  for 
the  ellipse.  The  pressure  on  the  windward  ray  of  the  elliptic-cone  waverlder 
Is  less  than  that  of  Its  circular-cone  counterpart,  but  the  pressure  at  the 
lips  Is  greater  owing  to  the  greater  Inclination  of'  the  shock.  There  are 
also  elliptic  cone  waverlders  with  zero  and  negative  dihedral,  corresponding 
to  the  circular-cone  waverlders  shown  In  Figs.  4.5  and  4.6,  but  these  will 
not  be  discussed  here. 

6.5  Further  Remarks 

Elliptic  waverlders  with  wedge-shock  upper  surfaces  can  also  be  constructed, 
such  as  was  done  for  the  circular  cone  configuration  shown  In  Fig.  4.7.  Because 
there  Is  an  additional  parameter,  e/s,  Involved  In  fitting  the  wedge-shock 
to  the  waverlder,  the  analysis  Is  very  involved  and  will  not  be  considered 
here. 

Other  characteristics  of  Interest  for  the  families  of  waverlders  presented 
here  are  lift,  drag,  and  stability.  The  pressure  coefficient  on  the  conical 
stream  surface  can  be  obtained  from  the  equation 

c Cp  (9)  Cp  (0)  . S (8) 

_£  (0,5)  * — 9 + (3.)  — 1 cos  + (1)  — i cos  2<|>  , (6.16) 

62  62  5 6 15  5 

which  is  In  hypersonic  similarity  form  and  can  be  regarded  as  known  [13,14], 

The  lift  and  drag  characteristics  can  thus  be  determined  without  great  difficulty. 

The  pitching  moment  can  also  be  determined,  but  general  stability  characteristics 

require  a knowledge  of  off-design  flow  conditions.  These  considerations  are 

fruitful  subjects  for  further  research. 
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7.  Some  Considerations  Regarding  Control  Surfaces 
7 . 1 Control  Flaps 

For  controlled  flight,  control  surfaces  are  necessary  for  changing  and 
maintaining  pitching,  yawing,  and  rolling  motions  as  well  as  aerodynamic  trim. 

The  waverlder  configurations  that  have  been  derived  so  far  generally  have 
horizontal  or  nearly  horizontal  surfaces  where  flaps  can  be  placed  to  produce 
nearly  vertical  forces  by  their  deflections.  If  enough  dihedral  is  present, 

It  may  be  possible  that  the  flaps  could  produce  the  behavior  of  a vee  tall 
and  thus  obviate  the  need  of  a vertical  fin.  In  general,  however,  the  design 
of  controlled  flight  necessitates  the  consideration  of  vertical  control  sur- 
faces. In  this  section  we  shall  formulate  the  analysis  of  vertical  fins  that 
are  themselves  waverlder  configurations. 

7 • 2 Vertical-Fin  Control  Surfaces 

We  consider  a wedge-shock  caret  waverlder,  such  as  shown  In  Figs.  2a  and 
2b,  but  that  Is  nonsymmetrlc,  that  is,  the  polar  angles  describing  the  free- 
stream  surfaces  are  unequal.  A typical  configuration  that  shall  be  utilized 
as  a vertical  fin  is  shown  in  Fig.  7.1.  Let  one  freestream  surface  be  denoted 
by  tan  Gg  = o$  (for  small  angles)  to  represent  the  angle  between  a freestream 

oo  oo 

cone  axis  and  the  corresponding  cone  shock  The  other  freestream  surface  which 
shall  represent  the  left  half  of  a vertical  fin  is  denoted  by  the  conical  angle 
<3y.  The  flow  deflection  angle  Is  denoted  by  A and  the  shock  angle  Is  denoted 
by  8W<  All  these  polar  angles  are  regarded  as  small  and  are  shown  as  projections 
on  a fixed  Z plane.  The  dihedral  angle  between  the  two  freestream  surfaces 
Is  denoted  by  ^ . 

From  geometrical  considerations,  the  law  of  cosines  yields 
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0S2  - 20s  0f  COS  + q\  = [(0S2  - e*)1A  + (of  - ej)!/2]2  . (7.1) 


Expanding  this  equation  and  solving  for  yields 


eS  °S  2 1/2 

of  ■ cos  '}„  t sin  t.  [(  J-Z-  ) - 1]  • (7-2) 


Since  ©s  /e^  Is  always  positive  the  plus  sign  must  be  chosen  when  ^ > 90° 

CD 

since  cos  >p  is  then  negative.  For  a continuous  variation,  a change  in  sign 

03 

occurs  when  the  radical  vanishes.  This  does  not  occur  since  ©s  /6W  > 1. 

oo 

Hence  the  plus  sign  Is  to  be  selected  in  Eq.  (7.2). 

Equation  (7.2)  follows  from  the  geometry  of  the  configuration.  The 
shock  conditions  of  gas  dynamics,  however,  require  that  £w  is  determined  by  the 
flow  deflection  a and  the  freestream  Mach  number  as  given  by  Eq.  (2.1).  We 
can  thus  rewrite  Eq.  (7.2)  as 


0^~  a cos  ^ + sin  [( 


2 


J 

B 


and  determine  <5/ 0W  from 


(7.3) 


K6  <7>‘ 


1 ] 


(7.4) 


Thus  (of/6)  is  a function  of  oc  /s,  <6>  y,  cos  ifr,  and  \/&,  where  6 is  taken 

JO 

as  the  characteristic  angle  of  the  cone  flow.  When  a = 0,  6, . takes  the 
Mach  angle  as  its  minimum  value,  and  takes  a minimum  value  given  by 


bf  . 

mm 


3 COS  ik  + sin  <(/  [(  — ) K:  - 1] 


i/: 


(7.5) 
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7 , 3 Vertical  Fin  for  the  El  1 Iptlc-Cone  Have rider 

We  consider  the  elliptic-cone  waverlders  shown  in  Fig.  5.4  with  the  con- 
ditions K.  = 1.3,  y = 1-4,  and  e/5  * 0.1.  In  this  case  we  match  the  freestream 

0 

upper  surface;  of  the  cone  waverider  with  the  freestream  surface  of  the  half- 
fin  analyzed  in  Section  7.2.  Wc  thus  set  = $s  and 

0$  3 0*  a b - eg  cos  2<|>s  . (7.6) 

SO 

For  the  lip  at  <es  * 100°,  we  get  for  various  values  of  A/6 


A/5 

0f/ 6 

0 

1.06 

0.1 

1.21 

0.2 

1.40 

Since  the  freestream  surfaces  are  matched,  the  vertical  fin  can  start  at 
the  cone  vertex,  or  at  some  other  position.  The  configuration  for  which  the 
fin  starts  at  the  vertex  is  shown  in  Fig.  7.2  for  a/6  ■ 0.2.  In  this  case  the 
fin  shock  is  attached  to  the  conical  lip.  The  case  where  the  fin  begins  half-way 
back  on  the  elliptic-cone  undersurface  Is  shown  In  Fig.  7.3  for  a/6  ■ 0.2.  In 
a real  flow,  owing  to  viscous  boundary-layer  shock  Interaction,  the  flow  would 
not  be  as  well-behaved  as  Illustrated  here. 

7.4  Vertical  Fin  for  Inclined-Cone  Waveri ders 

When  the  basic  conical  waverider  surfaces  are  derived  from  cones  at  angle- 
of-attack,  the  freestream  surfaces  are  inclined  at  an  angle  a with  the  geometric 
cone  axis.  The  formulas  in  Section  7.2  must  be  adapted  to  these  situations 
by  the  appropriate  evaluation  of  ^ and  9S  . The  unit  normal  vector  to  the  free- 

on 

stream  surface  passing  through  lip  shock  and  the  freestream  cone  axis  is  given 
by  hOT,  defined  by  Eq.  (6.14).  If  e denotes  the  unit  normal  vector  to  the 


vertical  plane,  then  we  obtain  * from  the  condition  cos  'L 


[a-  (f)(1-y)  cos  <p  - (£•)  g cos  2*  ] cos  * + f sin2$ 


J l>  - (|)(l-g)  cos  $s  - (J)  g cos  2$s]2  + (J)2  sin2  ^ 

When  a • 0,  this  reduces  to  the  result  <|>„  ■ it  - <ps . There  is  no  discrepancy 
in  this  result  compared  to  Section  7.3  If  It  is  realized  that  the  windward  ray 
was  taken  as  $ ■ 0 for  the  inclined-cone  waveriders,  whereas  the  non-inclined 
elliptic-cone  waverlder  utilized  $s  > it/2  for  the  conical  surface. 

From  geometrical  considerations,  we  also  obtain 

*2  * V* 

0$  a [0$  + - 20s  3 COS  $s]  , ( 

00 

where 


0S  ■ S + ag  COS  1>S  - 2eg  cos  2$, 


We  shall  not  consider  further  calculations  here. 


8.  Concluding  Remarks 


By  means  of  stream  surfaces  obtained  from  angle-of-attack  and  cross-section 
eccentricity  perturbations  of  the  basic  supersonic  axlsymmetric  flow  past  a 
circular  cone,  aerodynamic  lifting-body  configurations  have  been  derived.  The 
emphasis  has  been  on  a systematic  parametric  study  on  the  various  configurations 
that  can  be  obtained.  The  configurations  have  attached  shocks  on  sharp  leading 
edges  and  thus  can  be  described  as  conical  waveriders.  Utilization  of  wedge- 
shock  caret-waverlder  results  leads  to  a combination  of  configurations  and  to 
vertical-fin  control  surfaces.  Generalization  of  the  perturbed  flow  past  a 
circular  cone  to  allow  for  generalized  non-circular  cross  sections  [19]  can 
lead  to  further  aerodynamic  combinations. 

The  analysis  has  proceeded  within  the  framework  of  hypersonic  small  - 
disturbance  theory,  and  approximate  analytic  formulas  have  been  derived  that 
apply  over  a wide  range  of  conditions.  Although  pressure  distributions  have 
not  been  calculated,  the  pertinent  formulas  have  been  presented.  Further 
calculations  by  Interested  Investigators  can  be  performed  readily. 

The  lifting-body  configurations  that  have  been  presented  appear  attractive 
in  terms  of  high  lift-drag  ratio  requirements.  Further  work  Is  required  to 
account  for  other  aerodynamic  factors.  Some  of  these  are 

1)  lift,  drag,  and  moment  characteri sties , 

2)  boundary-layer  growth  and  related  viscous  effects, 

3)  off-design  effects, 

4)  details  of  flaps  and  other  control  surfaces, 

5)  unsteady  flow  and  dynamic  forces  and  moments, 

6)  blunted  edges  and  noses, 

7)  experimental  results. 
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These  subjects  are  part  of  a continuing  program  of  research  at  the  University 
of  Oklahoma,  sponsored  by  the  United  States  Air  Force. 
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APPENDIX  A 


FORMULAS  FOR  THE  BASIC  AXISYMMETRIC  CONE  FLOW 


Formulas  for  the  hypersonic  flow  field  In  the  shock  layer  of  an  unyawed 
cone  can  be  obtained  from  the  analysis  of  Rasmussen  [17].  They  are  also  derived 
alternatively  in  the  Appendix  of  reference  [11].  For  small  angles  within  the 
framework  of  hypersonic  small -disturbance  theory,  the  radial  and  polar  velocity 
components  are  given  approximately  by 
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and  s M^fi  Is  the  hypersonic  small -disturbance  parameter. 

The  temperature,  speed  of  sound,  and  Mach  number  are  given  In  terms  of 
the  velocity  for  homentropic  and  homenergic  flow.  We  obtain 
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where  errors  of  order  <52  have  been  ignored.  Since  the  flow  Is  homentropic, 
the  pressure  Is  determined  by 
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The  temperature  change  (and  hence  the  pressure  change)  across  the  shock  layer 


is  small,  and  hence  we  can  use  the  first  term  of  a binomial  expansion  to  write 
P„0)  „ T M - T„(a) 
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When  the  pressure  jump  across  the  shock  is  taken  into  account. 
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we  can  finally  write  the  pressure  distribution  across  the  shock  layer  as 


po(0)  , . r , , (y+1)  Ki2+  2 , , S2  , , e2  ,, 

"p — 3 1 + T C 1 TT-TF' Kj  +' 2”  { 1 • p-  + ln  |r  >1 

t»  Q 


(A. 8) 


*»-  • ••I’ 


* *«**•• 


72 


APPENDIX  B 

FORMULAS  FOR  THE  INCLINED  ELLIPTIC  CONE  FLOW  PERTURBATIONS 


The  flow  perturbations  outside  the  vortical  layer  are  a linear  combination 
of  the  angle-of-attack  contribution  and  the  cross-section  eccentricity  con- 
tribution, as  given  by  Eqs.  (3.6).  At  the  shock,  the  angle-of-attack  velocity 
perturbations  are  given  by  [11,12] 
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is  the  density  ratio,  and  a = s/>5  . The  shock  eccentricity  factor  is 


given  by 
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At  the  shock,  the  cross-section  eccentricity  velocity  perturbations  are 
given  by  [13,14] 
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and  the  shock  eccentricity  factor  is  given  by 
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The  pressure  disturbance  In  the  shock  layer  is  given  by  [11,12,13,14] 
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The  azimuthal  velocity  perturbation  w^a)  is  related  to  the  radial 
perturbation  u^e)  by  [13,14] 
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and  where  n r 1 for  the  angle-of-attack  perturbation  and  n = 2 for  the  cross- 
section  eccentricity  perturbation.  The  pressure  at  the  body  surface  0 * 6+0(a,e) 
does  not  depend  on  vL  since  vQ(s)  = 0.  Thus,  replacing  u 1 ( <s ) by  Wj ( 6)  by 
means  of  Eq.  (B.12)  yields  the  body  pressure  perturbation  as  a function  of  the 
azimuthal  velocity  perturbation: 
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When  Hq(5)  is  evaluated  by  means  of  Eqs.  (A.l)  and  (A. 2)  and  it  is  noted  that 
aQ2  (0)  varies  only  a small  amount  across  the  shock  layer,  such  that  aQ2  (0)=  aQ2  (e), 
then  H ( 6 ) has  approximately  the  value  [13,14] 
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Thus,  since  u„U)/V  ■“  1 + 0(<52)  , we  can  approximate  Eq.  (B.14)  by 
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The  coefficient  of  Fj  is  small  enough  to  neglect,  and  we  therefore  obtain 
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The  pressure  perturbation  at  the  surrace  is  thus  very  nearly  proportional  to  the 
azimuthal  veloi  ,ty.  This  result  is  consistent  with  the  fact  that  both  the 
pressure  and  first-order  azimutnai  velocity  are  uniformly  valid  across  the 
vortical  layer. 
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The  azimuthal  velocity  perturbations  evaluated  at  the  surface  have  the 
values  [11 ,12,13,14] 
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The  perturbation  pressure  variation  across  the  shock  layer  can  be  obtained 
from  Eqs.  (B.10)  when  u1 , ut,  Vj,  and  V]  are  evaluated  by  the  formulas 
of  references  [11,12,13,14].  Since  the  values  of  the  pressure  at  the  shock 
and  at  the  body  are  given  here,  a simple  linear  variation  with  0 across  the 
shock  layer  is  adequate  for  design  calculations. 
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